Time and frequency transfer with a microwave link 
in the ACES/PHARAO mission 



P. Delva*, F. Meynadier, P. Wolf, C. Le Poncin-Lafitte and P. Laurent 
LNE-SYRTE, Observatoire de Paris, CNRS et UPMC, 
61 avenue de 1' Observatoire, 75014, Paris, France 
*Email: Pacome.Delva@obspm.fr 



Abstract — The Atomic Clocks Ensemble in Space 
(ACES/PHARAO mission), which will be installed on board 
the International Space Station (ISS), uses a dedicated two-way 
Micro- Wave Link (MWL) in order to compare the timescale 
generated on board with those provided by many ground 
stations disseminated on the Earth. Phase accuracy and stability 
of this long range link will have a key role in the success of 
the ACES/PHARAO experiment. SYRTE laboratory is heavily 
involved in the design and development of the data processing 
software : from theoretical modelling and numerical simulations 
to the development of a software prototype. Our team is working 
on a wide range of problems that need to be solved in order 
to achieve high accuracy in (almost) real time. In this article 
we present some key aspects of the measurement, as well as 
current status of the software's development. 

I. Introduction 

The ACES/PHARAO mission is an international metro- 
logical space mission aiming at realizing a time scale of 
high stability and accuracy on board the International Space 
Station (ISS). Relative frequency stability (ADEV) should be 
better than a y — 10 -13 -r -1 / 2 , which corresponds to 3- 10 -16 
after one day of integration (see fig[T]); time deviation (TDEV) 
should be better than 2.1 • 10 -14 • r 1 / 2 , which corresponds 
to 12 ps after one day of integration (see fig(2]). Absolute 
frequency accuracy should be around 10 -16 . 

This mission is an international cooperation of more than 
150 people. PI laboratories are SYRTE/Paris Observatory, 
LKB/ENS and Neuchatel Observatory, and leading space 
agencies are the European Space Agency and CNES, the 
French space agency. Many industries are involved, the main 



ones being EADS/Astrium, TimeTech and Thales. All are 
working together to meet the scientific objectives of the 
mission: 

• Demonstrate the high permformance of the atomic clocks 
ensemble in the space environment and the ability to 
achieve high stability on space-ground time and fre- 
quency transfer. 

• Compare ground clocks at high resolution on a world- 
wide basis using a link in the microwave domain. In 
common view mode, the link stability should reach 
around 0.3 ps after 300 s of integration; in non-common 
view mode, it should reach a stability of around 7 ps after 
1 day of integration (see fig(2]). 

• Perform equivalence principle tests. It will be possible 
to test Local Lorentz Invariance and Local Position 
Invariance to unprecedent accuracy by doing three types 
of tests: a test of gravitational red- shift, drift of the fine 
structure constant and of anisotropy of light. 

Besides these primary objectives, several secondary objec- 
tives can be found in [1]. For example, if the theory of general 
relativity is considered as exact, then the measurement of 
gravitational redshifts can be used to measure gravitational 
potential differences between different clock locations. It is 
a new type of geodetic measurements using clocks called 
relativistic geodesy. 

In this article we describe in details the Micro- Wave Link 
(MWL) used in the ACES/PHARAO mission, and developed 
by TimeTech (TT). First we describe one-way and two- 




Fig. 1. PHARAO (Cesium clock) and SHM (hydrogen maser) expected 
performances in Allan deviation. 



Fig. 2. Performance objective of the ACES clocks and the ACES space- 
ground time and frequency transfer expressed in time deviation. 



way links theoretically and introduce the SYRTE Team (ST) 
observables. In the second part we describe how works the TT 
modem, and what is the link between the modem observables 
(TT observables) and the ST observables. Finally we present 
the status of the data analysis software and of the simulation 
we are developing. 
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II. The Micro-Wave Link (MWL) 

The Micro- Wave Link (MWL) will be used for space- 
ground time and frequency transfer. A time transfer is the abil- 
ity to synchronize distant clocks, i.e. determine the difference 
of their displayed time for a given coordinate time. The choice 
of time coordinate defines the notion of simultaneity, which 
is only conventional. A frequency transfer is the ability to 
syntonize distant clocks, i.e. determine the difference of clock 
frequencies for a given coordinate time. Here we suppose that 
all clocks are perfect, i.e. their displayed time is exactly their 
proper time. Proper time r is given in a metric theory of 
gravity by relation: 



c 2 dr 2 



-g ai3 dx a dx (3 , 
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where g a p is the metric, c the velocity of light, {x a } the 
coordinates and Einstein summation rule is used. We use in 
this article the notation [.], which is the coordinate / proper 
time transformation obtained from eq.Q, and = tj — ti 
for coordinate time interval^] 

The MWL is composed of three signals of different fre- 
quencies: one uplink at frequency ~ 13.5 GHz, and two 
downlinks at ~ 14.7 GHz and 2.2 GHz. Measurements are 
done on the carrier itself and on a code which modulates 
the carrier. The link is asynchronous: a configuration can be 
chosen by interpolating observables. In the following we give 
a formal description of one-way and two-way links for code 
observables. The principle for carrier observables is the same 
except that periods cannot be identified, leading to a phase 
ambiguity. 

A. One-way link 

1) Experiment: let's consider a one-way link between a 
ground and a space clock represented respectively by subscript 
g and s. The sequence of events is illustrated on fig[3] At 
time coordinate ti, clock g displays time T\ and modem Mg 
produces a code C 1 . This code modulates a sinusoidal signal 
of frequency / and sent at coordinate time £2 by antenna g. 
The delay between the code production and its transmission 
by antenna g is A 9 = [T12] 9 , expressed in local frame of 
clock g. Antenna s receives signal C 1 at coordinate time £3, 
and transmit it to modem Ms and clock s which receives it at 
coordinate time £4, with a delay A s = ^34]*. Clock s displays 
time ri and modem Ms produces the code C 1 at coordinate 
time £5. 

^.g. [T12] A is the transformation of coordinate time interval T12 in proper 
time of clock A, and [Ar A ]* is the transformation of proper time interval of 
clock A in coordinate time t. 




(b) Space-time diagram 
Fig. 3. Schematic representation of the one-way link. 



We use superscript g or s on proper times r for clocks g 
or s, and we express proper time as a function of coordinate 
time. Then we can write 



Tl = T»(*l) = T*(* 5 ). 



(2) 



We define the ST (SYRTE Team) observable Ar s given by 
modem Ms with: 



Ar s (r s (U))=r s (t 5 )-r s (U). 



(3) 



This observable is dated with proper time of clock s when this 
clock receives code C 1 from antenna s. It can be interpreted 
as the difference between the time of production of code C 1 
by clock s, and time of reception of same code C 1 sent by 
clock g, all expressed in proper time of clock s. 

2) Desynchronisation: desynchronisation between clock g 
and s is written in an hypersurface characterized by coordinate 
time t = constant. From eqs.([2])-([3} it is straightforward to 
deduce it for coordinate time t±\ 



t 8 (U)-t°(U) = -Ar s (r s (U))- T 23 + [A 9 + A s ] 



(4) 



Similar formulas can be obtained for the desynchronization 
at coordinate times t\ and £5. This expression has been 
obtained for the uplink, from ground to space. To obtain the 
desynchronisation with downlink observables, all you have to 
do is replace g and s in eq.Q. 

B. Two-way link 

1) Experiment: let's consider now a link which is com- 
posed of two one-way links, between a ground and a space 
clock represented respectively by subscript g and s. The two 
sequences of events are illustrated on fig. [4] The uplink (from 
ground to space) has a frequency j\ and is repre sented by 
coordinate time sequence 

(t?,ti,t2, *§,*?) (figpOOl- This link 
is defined with the relation: 

^(t?)=T S (t?). 
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(a) Uplink sequence of events 
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(b) Downlink sequence of events 




(c) Space-time diagram 
Fig. 4. Schematic representation of the two-way link. 



The downlink has a frequency f 2 and is repre sented by 
coordinate time sequence (£3, £3, £4, £4, £g) (fig 3(b) ). This link 
is defined with the relation: 



2) De synchronisation in a two-way configuration: a two- 
way configuration is defined by r 9 {t\) = r s (t^), i.e. the code 
C 1 of link fi is the same as code C 2 of link f 2 , and they are 
locally produced and sent at the same time: t\ = £g (at clock g) 
and £3 = £7 (at clock s). Then we calculate desynchronisation 
between clocks g and s at coordinate time t\ as: 



1 r 
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Arf (r^°))+Af + A| 
At s (T s (t° 2 )) + Af + Af. 



(5) 



(6) 



3) Desynchronisation in a A configuration: in the 
ACES/PHARAO mission we use the so-called A configura- 
tion. This configuration minimizes the error coming from the 
uncertainty on ISS orbitography (in [2] it has been shown that 
in this configuration the requirement on ISS orbitography is 
around 10 m). The A configuration is defined by £2 = £3, i.e. 
code C 2 is sent at antenna s when code C 1 is received at 
this antenna. This configuration is obtained by interpolating 
the observables. Then it can be shown that desynchronisation 



between clocks g and s at coordinate time £2 is: 

r«(t 2 ) - r»(«2) = \ (Arl(^) - A< (t°) + [T 34 - T 12 f) . 

(7) 

C. Approximations 

1 ) Coordinate / proper time transformation: in equation ([7]) 
remains one transformation from coordinate to proper time. 
We know that T12 ~ T34 ~ 1 ms. During this time interval 
we can consider that the gravitational potential and velocity 
of the ground station are constant. Therefore we can do the 
approximation: 



[T 34 - T 12 } 9 = (1 - e g {t 2 )) (T 34 - T 12 ) , (8) 



where 



GM ^(£) 



r^(£)c 2 2c 2 ' 



M is the Earth mass, r g (t) and v p (£) are the radial coordinate 
and the coordinate velocity of the ground clock at coordinate 
time £. Orders of magnitude of these corrective terms are: 



GM r 
r g c 2 

*,2 



34 



0.6 ps 



^T 34 - 0.002 ps 

The gravitational term is just at the limit of the required 
accuracy. The velocity term is well below, so we can neglect 
it. Final formula for desynchronisation is then: 

1 



r s fe)-r^(t 2 )= 2 (Ar^)-Ar^) 



1 



GM 



(T 34 -T 12 ) . 



(9) 



2) A configuration: in the A configuration we suppose that 
T23 = 0. However, this will never be exactly and it will be 
known with a precision £T23. This will add a supplementary 
delay S(r s — r 9 ) to desynchronisation ([7]): 

^(r s -^)(£2) = (e,(£2)-e s (£ 2 ))^T23 
Orders of magnitude are: 
GM ( 1 1 

2 2 
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With the required accuracy on the MWL, \S(r s —r 9 ) \ < 0.3 ps, 
we deduce the following constraint on ^23: 

^23 < 0.9 ms. 

This constraint is much less constraining than the one coming 
from orbitography, which is ST23 < 1 /is (see [2]). 



D. Atmospheric delays 

The downlink is composed of two one-way links of fre- 
quencies /2 and fs, represented respectively by coordinate 
time sequence (£3, £3, £4, £4, £g) and (£5, £5, ^5 tg). These 
two links are affected by a ionospheric delay that depends on 
their respective frequencies, whereas the tropospheric delay 
does not depend on the link frequency (we neglect dispersive 
effects). We write: 

m ^34 
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where R^ = \xR X (tj) — xtx(U)\ is me range, xt x and xr x are 
respectively position vectors of space and ground antennas, 
rj x \x Tx \ and r Rx |f Rx |. 

Ionospheric and tropospheric delays are around or below 
100 ns, whereas Shapiro delay (term in c -3 ) is below 10 ps 
for the ACES/PHARAO mission (see [3] and figjS]). 

1) Ionospheric delay: in order to deduce ionospheric de- 
lays, we combine the two ground observables to be free of 
tropospheric delays. We obtain: 

Ar'(r'(tg)) - AT°(T°(ft)) = [T 34 - T 56 ] s 
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(12) 



Here we impose that T® 6 = 0, ie. both signals are sent by 
clock s at the same time. However, this will never be exactly 
zero, there will be a remaining ST® 6 introducing a timing error 
ST ~ (e s (tl) - e g (f A )) JT 4 ° 6 . With a required accuracy ST < 
0.3 ps, we obtain the following constraint: 

8T% 6 < 0.9 ms. 

We expect that |T 34 - T 56 | < 100 ns (see [3]); therefore 
we can neglect the coordinate to proper time transformation 
in eq.([T2|). We can also neglect this transformation for the 
delays. Then eq.([T2]) is equivalent to: 

ArKO - Ar£ (i°) = T 34 - T 56 (13) 

From eqs.([T0|), fTT] ) and ([13) we obtain: 

R34 — R56 



*56 



'(/ 3 ) 



^34 



\f 2 ) = Ar£ (i°) 



(14) 

where we neglected the difference of the Shapiro delays 
between the two downlinks, which can be shown to be 
completely negligable. 

Now we can calculate S the Slant Total Electron Content 
(STEC). The ionospheric delay affects oppositely code and 
carrier and may be approximated as follows: 
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(15) 



B-k)dL (16) 



where N e is the local electron density along the path, STEC 
S = J N e dL, B is the Earth's magnetic field and k the unit 
vector along the direction of signal propagation. It has been 
shown that higher order frequencies effect can be neglected 
for the determination of desynchronisation [3]. 

We suppose that for a triplet of observables 
{At 3 (t 3 (4)), Ar^r^(^)), Ar^r^tg))}, variations of 
the direction of signal propagation and of magnetic field 
along the line of sight does not change: B ~ Bq. Then: 
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where #0 is the angle between B and the direction of propa- 
gation of signal and f%. Then we obtain: 
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[A^(/ 3 )-A 3 T(/2)] ca = 
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(20) 

These equations, together with equation ([14]), give the 
STEC S. The value of S can then be used to correct the 
uplink ionospheric delay. 

2) Tropospheric delay and range: by adding ground and 
space observables of links fi and ji we obtain: 



At s (t s (^)) + At2(t^)) + A 



[Af + Aif 
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As in the previous section, it can be shown that [T^Y — 
[ T 2sY = if T 2°3 is known with a precision ST$ 3 < 0.9 ms. 
We neglect the coordinate to proper time transformations for 
delays and obtain: 
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Then, neglecting Shapiro time delays we obtain: 

R12 + R3A (. GM 



- (A™°(/i) + A™°(/ 2 ) + A7° + A*° po ) 

This equation shows that range and tropospheric delays are 
degenerated. Range can be calculated with a model for tropo- 
spheric delay, and tropospheric delay can be calculated from 
an estimation of range. 
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Fig. 5. Both emitter's and receiver's signals are represented against coordinate time scale. Red dots on the filtered beatnote indicate zero crossings on 
ascending edge. n m is the number of red dots between two 12.5 PPS pulses. T m _i is the proper time of the first red dot in the same sequence. Propagation 
times are represented by black arrows. Here signal frequency is much lower than in reality, and Doppler effect is strongly magnified in order to show variation 
of n m . 



III. Micro-Wave Link modems 

We explain here basic principles of the ground/space 
modems developed by TimeTech/Astrium for the 
ACES/PHARAO mission, that will be linked to the clocks 
and the antennas. This principle is illustrated on fig (5] At 
emitter and receiver are generated a PPS signal (one Pulse 
Per Second), a 12.5 PPS (one pulse every 80 ms, the period 
of measurements), and a periodic signal (either code at 
100 MHz or carrier). Let e be the emitter and r the receiver. 
A PPS signal sent at local time r p e ps of the emitter is received 

at local time of the receiver. Local time is recorded 

pps 

by the modem for each received PPS. 

When received, the periodic signal (blue) is mixed with a 
local oscillator (yellow) which frequency is not far from the 
received frequency, and filtered to obtain the low frequency 
part of the beatnote (green). The beatnote frequency is around 
195 kHz for code and 729 kHz for carrier. The receiver 
modem records the time of the first ascending zero-phase of 
the beatnote signal after the 12.5 PPS signal. We call this 
observable T m , where m is the number of the 80 ms sequence. 
Finally, the modem counts the number of ascending zero-phase 
n m during sequence m. 

Tp ps , T m , rim and m are the basic observables of the modem, 
called TT observables, and are recorded for code and carrier 
signals. The modem internal clock is reset every 4 s. However, 
code observables can be linked to UTC time, which permits 
to solve the phase ambiguity between each passage for code 
observables. 



A. From TT to ST observables 

Let (j) e (i~ e ) and r (r r ) be respectively the phase of emitted 
and received signals, changing with local time of emitter and 
receiver. Let's consider two signals: one emitted at emitter 
local time Tq and received at receiver local time Tq, and 
another one, emitted at rf and received at t[ (see fig. [6]). The 
phase increase between these two signals is equal at emitter 



(22) 



and receiver: 

^e(rf) - </>eK) = <Mt[) - MrS). (21) 

The received signal is mixed with a local oscillator signal such 
that the beatnote phase is: 

r ) -(/V(r r ) (for code) 
^L.o.(^ r ) (for carrier) 

Signs are different for code and carrier; we will write subse- 
quent formulas in a compact way with the sign in black for 
code and red for carrier. We assume <p(r) = cjr + cst, where u 
is the pulsation of the considered signal. Then, from eqs.([2T]) 
and ( [22] ) we deduce: 

+ 1 



Mr r ) 



J 0L.o.(r r 



VL.O. 



-(MtD-MtS)) (23) 



We introduce the ST observable, which links local time of 
emission to local time of reception: Ar r (r r ) = r e — r r . From 
eq.([23]) we get: 



Arf(Tf) - Ar£(r r ) = 



WL.O. 



1J (r[ - tS) 
(Mr[) ~ <h(r r )) 



(24) 
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Fig. 6. Link between phase, emitter and receiver local times. 



Let's apply this formula to the MWL modem by introducing 



the TT obervables: 



T m _i and rf 



From 



definition of T m and n m observables we know that 0b (T m ) 
0b(T m _i) = 27rn m . Then we deduce from eq.([24|): 

Ar m (T m ) = Ar^ l _ 1 (T m _i) 



+ 27rn m 



1 (T* 



l) 



(25) 



UJ e 



where Ar^ is the ST observable corresponding to se- 
quence m. This recursive formula allows to find all ST 
observables from TT observables, if the first term Atq (T ) 
is known. We notice than in case of zero Doppler, the last 
two lines of the equation cancel, and the ST observable Ar r 
should be constant with local time. 

Relative accuracy of ST observables during one passage of 
ISS is: 



195 kHz 
100 MHz 



729 kHz 



10 ns ~ 20 ps (code) 



v 13 5 GHz 10 n s ~ 0.5 ps (carrier) 

where accuracy of T m observables, ST m ~ 10 ns, can be de- 
duced from modem internal clock, which frequency is around 
100 MHz. However, ST m is underestimated here because other 
noise sources than the internal clock may count. The goal here 
is not to do a precise accuracy budget but rather get a lower 
limit. 

B. Initial term determination 

What is the first term of the iterative series ([25])? Our 
goal is to determine Atq(To) with an absolute accuracy 
< 100 ps, which is required for ground-space time transfer. 
For frequency transfer, we should be able to bridge the gap 
between two passages with an accuracy depending on the 
duration between them, which can be read on fig|2](e.g. 1.5 ps 
for two passages separated by one orbital period). 

Let 



Ar r (r r ) = r 

' pps V ' pps / ' i 
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(26) 



be the ST observables linked to TT observables r pps , and r[ ni 
be the receiver local time of generation of the PPS, which 
is by definition of the experiment equal to the local time of 
emission of the same PPS signal r pps . As r pps can be guessed 
from a UTC tag, Ar^ s (r^ s ) is known and can be linked to 
Ar r (r r ) with the help of eq.(24): 



Arff(rff)= Ar; s (r; ps ) + 



UJ e 



'ppsy 



--(^(r;)-$ b (r; s )) 

(27) 

From the values of T m and n m one can determine the 
beatnote phase at T pps . The uncertainty in that determination is 



5$b(Tp ps ) ~ co>5 • STm ~ 0.012 rad. The uncertainty of the two 
last terms on the rigth side of the equation is then ~ 20 ps, 
which is sufficient. However, ArJL is known with the modem 

ppr> 

internal clock accuracy, i.e. r5(AT pps ) ~ 10 ns. Even averaging 
on a complete data set is not sufficient, reaching ~ 580 ps 
accuracy with 300 points. Then we need a method to obtain 
a precise PPS observable. 

1 ) Precise PPS observable: the emitter PPS is phase coher- 
ent with the code phase and the receiver PPS is phase coherent 
with the local oscillator phase. Then we can use the internal 
counter and the code phase observables to follow the phase 
precisely and derive a more accurate pps observable that we 
call AT p r pp S . 

The receiver local oscillator is phase coherent with the 
receiver PPS (it has a zero crossing at rf ni ), then: 



$ L .0.(^) = 27TiV L .o. 



(28) 



where 27V L . . is an integer. Similarly the emitted (and received) 
code is phase coherent with the emitted (received) PPS, so we 
have 



* r (0 = 2ttN b 



(29) 



'pps/ -rvppsy 

where N B is an integer. Using eqs.([22l>, ([26]>, ( |28l > and <[29j 
we obtain: 

A^ s (t; s ) = ^ (2irN - Mt£ ps )) , (30) 

where 2N = 2(7Vl.o. — Nb) is an integer. Provided we can 
determine the integer 2 N exactly, the uncertainty on ATp pps 
is ^b/W.o. ~ 20 ps. The value of N is determined by using 
the direct measurement ATpp S (Tp ps ) in eq.pO]) above: 

1 



N 



(uj uo AT p r p S (Tp r p S ) + $b(r p r ps )) 



(31) 



Uncertainty of the second term in ([31} is ~ 0.012/(27r) ~ 
2 • 10 -3 , therefore negligible. However, uncertainty of the 
first term is ~ uj^o. • 10 ns/(27r) ~ 1, which is insufficient. 
One solution is to average over all PPS measurements of a 
continuous passage, which should be sufficient for realistically 
useful passages (e.g. > 50 s). 

Finally using the precise PPS observable AT ppps (r pps ) cal- 
culated from eq.([30|) in eq.([27|) reaches the required accuracy 
for time transfer. 

2) Bridging the gap: to perform ground-space frequency 
comparisons on long term one needs to link observables from 
two different ISS passages. It is then necessary to determine 
precisely the time elapsed from one passage to another thanks 
to the TT observables. Let's call Ato(to) and Ato(tq) the 
initial terms of two different passages (in this section we omit 
the r indices on ST observables Ar r and on local time r r ). 
The error on their absolute determination is ~ 20 ps. However 
we want to reach S (Atq(tq) — Ato(to)) < x, where x is a 
specification that depends on the duration separating the two 
passages and can be read from fig|2](e.g. 1.5 ps for one orbital 
period separation). From eq.p4|) we deduce: 

Mtq) ~ 0b(r o ) = -uJeiAr^) - Ar (r )) 

(32) 
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Fig. 7. Short term stability (TDEV) of ST observables, derived from either code/carrier phase or accumulated code/carrier phase. 



Moreover, we know that fa (tq) — fa(To) 
is an unknown integer. We deduce that: 



27rN g9 where N g 
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(33) 



(-^ e (AT^)- Aro(ro)) 

+ (^L.O. ~ UJ r ) (Tq - T )) 

It can be shown that the accuracies of the different terms in this 
equation are sufficient to determine N g whithout ambiguity. 



Then we deduce: 
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This equation links ST observables from two passages to 
~ 20 ps, which is not sufficient as can be seen from fig [2] 
Using several pairs of code observables to determine this 
quantity will not increase the accuracy, as the error for each 
data pairs will be correlated. Two solutions can be envisioned. 
One can use carrier phase observables. However, it remains to 
be seen if the phase ambiguity (integer N g ) can be solved for 
carrier phase. This will be studied in another article. Second 
solution would be to use another observable from the MWL 
modem: the accumulated phase T™°, which is the sum of all 
dates of ascending zero-phase during sequence m. 

In fig|7]is shown the short-term stability of ST observables, 
calculated either using T m or T^ c observables, for code 
(fig 7(a)[ ) and carrier (fig 7(b)). It can be seen that using 
observables increases measurements accuracy. However, code 
accuracy is still not sufficient to bridge the gap between 
two passages: we assume that ST observables uncertainty is 
5 (At) ~ 3 x TDEVo, where TDEVo is the Time deviation for 
r = 0.08 s, the period of measurements. Then, from fig. |7(a) 
5 (At) ~ 90 ps for code phase and 5 (At) ~ 18 ps for 
accumulated code phase. This is not sufficient to bridge any 
gap. From fig. |7(b) we deduce 5 (At) ~ 0.9 ps for carrier 
phase and 5 (At) ~ 0.15 ps for accumulated code phase. Any 
of these two observable is sufficient to bridge a gap of 30 mn or 
larger. It remains to be seen how to solve the phase ambiguity 
for carrier phase, which will be studied in another article. 

IV. Simulation and data analysis 



mission data. This software is written in Python language. In 
order to test it, we wrote a simulation that generates (noisy) 
TT observables, as well as theoretical ST observables. This 
simulation is written in Matlab language, and is as much as 
possible independant from the data analysis software. 

A. Simulation 

The simulation takes as input orbitography of the ISS and 
one Ground Station (GS) in a Celestial Reference System. 
From orbitography files it simulates the proper times given 
by ISS and GS clocks, and the time transfer between these 
two clocks, using a modelization of the MWL. Observables 
are given in terms of TT and ST observables. Moreover, 
theoretical values of the scientific products are given in order 
to test the data analysis software. 

In fig (8] we have plotted different contributions included in 
the signal time-of-flight of signal fi ~ 13.5 GHz. Minimum 
elevation of ISS is taken as 10 , and atmospheric parameters 
are temperature T = 298 K, pressure p = 1 bar and water 
vapor pressure e = 0.5 bar. A sinusoidal variation is added 
to atmospheric parameters, with a period of 24 hours, and a 
Chapman layer model is used to calculate the STEC. 
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The SYRTE team writes an independant data analysis 
software, in order to make the most of the ACES/PHARAO 



Fig. 8. Contribution of atmospheric and Shapiro delays in the simulated 
time-of-fligth of the fi signal. 
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Fig. 9. Pre-processing software: comparison between ST observables 
calculated from simulated TT observables, and theoretical ST observables, 
for carrier (left) and code (right). 



Tropospheric delay is dominant in the time-of-flight, with a 
value of several 10 ns. We used here a Saastamoinen model, 
which is not really reliable at low elevation of ISS. The 
dispersive part of troposphere has not been taken into account, 
and it remains to be seen if this is necessary. Ionosphere is 
dispersive such taht the ionospheric delay can be separated 
in two contributions: an effect that scales with 1/f 2 and one 
that scales with 1/f 3 (see eqs.([T5])-([T6|)). The second order 
contribution is around 1 ns and the third contribution around 
0.1 ps, below mission accuracy. However these effects are 
much larger for the = 2.25 GHz signal: around 40 ns for 
second order term and 20 ps for third order term. Then third 
order terms cannot be ignored. Finally the Shapiro delay of 
several ps is sligthly over the required accuracy. 

B. Data analysis 

An independant pre-processing software has been written, 



using equations from sec III It takes TT observables from the 
simulation, transform them to ST observables and compare 
the result to the theoretical ST observables coming from the 
simulation. One example can be seen on fig(9] It can be seen 
that ST observables are well recovered, with a noise which is 
coherent with previous estimations. However here the noise is 
underestimated because all noise sources has not be included 
(only the modem internal clock noise). The absolute value of 
the ST code observable is found thanks to the method of the 
initial term determination, with an uncertainty less than 20 ps 
(explaining why the data cloud is not centered on 0). The 
phase ambiguity for carrier observable has not been solved, 
explaining the constant bias between the recovered and the 
theoretical ST carrier observables. 

The full data analysis software is being written but not 



finished yet. We explain its basic principle on fig [10] A special 
care is taken for file naming, data classifying, file formats and 
conventions. Indeed many data from several different sources 
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Fig. 10. Full data analysis software: illustration of inputs and outputs. 



will have to be used and these issues can be critical. The 
software design has been build in a modular way, and now 
most of the building blocks are written. 

V. Conclusion 

We have written a theoretical description of one-way and 
two-way satellite time and frequency transfer and devel- 
oped a model of the Micro-Wave Link in the frame of the 
ACES/PHARAO mission. This description has been used to 
write a data analysis software and a simulation to test it. The 
simulation is written in its first version, and used to assess 
our pre-proceesing software. The design of the data analysis 
software has been done in a modular way, and most of the 
building blocks are ready. 

Several questions remains: how to solve the phase ambiguity 
for the carrier observable, what is the dispersive effect of the 
troposphere? 
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